AIAA Guidance, Navigation and Control Conference and Exhibit 
18-21 August 2008, Honolulu, Hawaii 


AIAA 2008-6504 


An Integrated Approach to Damage Accommodation in Flight Control^ 


Jovan D. BoskovicJ Nathan Knoebef and Raman K. Mehra'J 
Scientific Systems Company, Inc. 

Irene Gregory 11 

NASA Langley Research Center 


In this paper we present an integrated approach to in-flight damage accommodation in flight control. The 
approach is based on Multiple Models, Switching and Tuning (MMST), and consists of three steps: 

• In the first step the main objective is to acquire a realistic aircraft damage model. Modeling of in-flight 
damage is a highly complex problem since there is a large number of issues that need to be addressed. One of 
the most important one is that there is strong coupling between structural dynamics, aerodynamics, and flight 
control. These effects cannot be studied separately due to this coupling. 

• Once a realistic damage model is available, in the second step a large number of models corresponding to dif- 
ferent damage cases are generated. One possibility is to generate many linear models and interpolate between 
them to cover a large portion of the flight envelope. Once these models have been generated, we will implement 
a recently developed Model Set Reduction (MSR) technique. The technique is based on parameterizing damage 
in terms of uncertain parameters, and uses concepts from robust control theory to arrive at a small number 
of ’’centered” models such that the controllers corresponding to these models assure desired stability and ro- 
bustness properties over a subset in the parametric space. By devising a suitable model placement strategy, the 
entire parametric set is covered with a relatively small number of models and controllers. 

• The third step consists of designing a Multiple Models, Switching and Tuning (MMST) strategy for estimating 
the current operating regime (damage case) of the aircraft, and switching to the corresponding controller to 
achieve effective damage accommodation and the desired performance. 

In the paper present a comprehensive approach to damage accommodation using Model Set Design, MMST, 
and Variable Structure compensation for coupling nonlinearities. The approach was evaluated on a model of 
F/A-18 aircraft dynamics under control effector damage, augmented by nonlinear cross-coupling terms and 
a structural dynamics model. The proposed approach achieved excellent performance under severe damage 
effects. 


I. Introduction 

There has been an increased interest in addressing flight control problems of aerospace vehicles under adverse con- 
ditions, which are defined as failures, damage and/or upset. This interest is motivated by both the desire to ensure safe 
flight in the presence of adverse conditions and an increased terrorist threat from the man-portable air defense systems 
(ManPADS). There is a renewed interest in development of effective adaptive reconfigurable control systems that can 
recover from an adverse event such as damage and reclaim as much vehicle nominal performance as is physically pos- 
sible. The design of such control systems is a difficult and challenging problem since damage results in both structural 
and aerodynamic disturbances. For instance, when a wing is damaged, the resulting dynamic disturbance is primarily 
due to reduced wing stiffness and loss of resistance to flutter. The aircraft can continue to fly after damage only if either 
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the static stiffness is sufficient to withstand aerodynamic loads, or the damage does not cause loss of flutter control- 
lability and wing destruction. 16 In addition, wing damage generates asymmetric disturbances that couple longitudinal 
and lateral dynamics. Hence the general problem of aircraft damage modeling and accommodation is highly complex, 
and there are only a few results that address this problem under highly restrictive simplifying assumptions. 

The problem of controlling a flexible aircraft with significant structural dynamics is highly complex. Initial analy- 
ses were performed and solutions to this problem were proposed in the literature. 7 ' la 12-14 To address the problem of 
control effector damage, Wise et al 19 assumed that the control effector damage causes the changes in the gain effective- 
ness matrix, and generates state-dependent disturbances. The scheme estimates the gain effectiveness matrix using a 
standard adaptive control approach, and the unknown damage-generated nonlinearity using radial basis function neural 
networks. It was assumed that there is a single nonlinearity due to damage, and that the bounds on the neural network 
weights are known a priori, which results in the boundedness of the signals in the system. The latter assumption may 
not be realistic since prior knowledge of the effect of damage is not generally available. In the area of wing damage 
accommodation, Nguyen and Krishnakumar 18 developed a vortex-lattice based damage model for various damage con- 
figurations, and also added the effect of the loss of symmetry of the damaged aircraft. A control design model was 
developed through linearization, and a direct adaptive control technique was used to compensate for the damage. The 
model does not take into account structural model changes and was evaluated for relatively small amounts of damage. 
The approach by Boskovic and Mehra 6 is based on a restrictive assumption regarding the damage model. The main 
assumption is that different levels of damage can be modeled as matrix polynomials in the damage paramater that 
belong from zero to one. No structural effects of damage were taken into account. Similar assumptions were made by 
Boskovic, Bergstrom and Mehra, 2 * where the simultaneous effect of actuator failures and the control effector damage 
was considered. 

In this paper we present an efficient Integrated Damage & Adaptive Control System (IDACS) that accommodates 
control effector damage. A block diagram of IDACS is shown in Figure 1. The IDACS is based on the Multiple 



Figure 1. Structure of the Integrated Damage-Adaptive Control System (IDACS). 


Models, Switching and Tuning (MMST) technique that has already been applied to the problem of actuator failure 
accommodation. 5 ' 17 The current approach is based on developing coupled structural dynamics and aerodynamic mod- 
els for different damage amounts, deriving a suitable parametrization of the models in order to generate a relatively 
small number of models covering the entire parametric set, and building corresponding multiple-model estimators. The 
estimates of unknown nonlinearities can be generated using either neural networks or Variable Structure techniques. If 
a solution to the wing damage accommodation problem exists, the use of these estimates in the reconfigurable control 


2 of 25 


American Institute of Aeronautics and Astronautics 




law will assure stabilization of the aircraft in the presence of wing damage. This is an important problem that has not 
been solved so far, and its solution will contribute to improved flight safety of both commercial and military aircraft. 

This research has been funded by the Integrated Resilient Aircraft Control (IRAC) project under NASA’s Aviation 
Safety Program. The goal of the IRAC project is to arrive at a set of validated multidisciplinary integrated aircraft 
control design tools and techniques for enabling safe flight in the presence of adverse conditions (ex: faults, damage 
and/or upsets). Major elements of the IRAC project involve integrated modeling of aerodynamics and structures; 
adaptive control under adverse conditions; and verification and validation of models and control methodologies. Hence, 
the approach to mitigate the impact of damage on the stability and performance of aircraft proposed in this paper 
contributes to two of the three areas emphasized in IRAC in an integrated fashion. 

To conduct the research under this project we obtained a simple model of wing damage in F/A-18 aircraft form 
Boeing that simulates control effector damage effects on aircraft dynamics, and augmented it with cross-coupling 
terms due to assymmetry caused by damage, as well as by a structural model of the TAMU wing. 20 In the paper 
present a comprehensive approach to damage accommodation using Model Set Design, MMST, and Variable Structure 
compensation for coupling nonlinearities. Several issues have been specifically addressed: (i) How to arrive at a simple 
Model Set Design (MSD) procedure? (ii) How to include the cross coupling terms into post-damage aircraft dynamics, 
and how to compensate for these terms? and (iii) How to integrate rigid-body and structural dynamics into a meaningful 
model and design a corresponding MMST controller? 

II. Damage Models and Baseline Control Strategies 

In this section we describe three damage models of increasing complexity that were used under the project. As- 
sociated with each of the models is a baseline control strategy used as basis for adaptive and MMST-based control 
design. 

Under this project we obtained a simple linear simulation model of wing damage in F/A-18 aircraft from Boeing. 
The flight regime is 0.7 Mach and 20, 000, and the corresponding A and B matrices are modified to include perturbations 
due to damage. This model was used as a baseline, and was subsequently augmented by inter-axes coupling effects due 
to damage, and by a structural model consistent with that derived for the TAMU wing. 20 

The nominal F/A-18 linear model is of the form: 

x = A n x + B n u (1) 


with perturbed states and inputs 

x - [u, q, , 6 , a, h, f3, p, r, <p, i//] T 

u = [ LEF/, LEF,., TEF,, TEF r , AIL,, AIL r , STAB,, STAB r , RUD,, RUD,., PLA,, PLA,f. 

In the case with no damage, the baseline controller is of the form: 

u = B\(—A n x + A m x + B m r ), (2) 

where B), denotes a pseudo-inverse of B n . 

The damage model includes the nominal system matrices A„ and B n , plus two additive damage perturbations ma- 
trices such that 

x — ( A n + AA( 6 \))x + ( B n + AB( 62 ))u (3) 

where AA = [Oioxi $102 Oioxs], ai is the second column of A, and A B = Oio X 9 ], and b\ is the first column of B. 
Boeing engineers feel that this model adequately describes small wing damage in F/A-18 aircraft. 

We note that the assumed damage perturbations are defined using two uncertain terms, namely 9\ and 62 , where 
only the columns of the state and input matrices associated with q and LEF, are affected. For this example it is assumed 
that all states are measured, actuator dynamics are negligible, and the control objective is to follow the reference model: 

Xm — A m x m + B m r. (4) 

In this case, the baseline controller is of the form: 

u = (B„ + AB 02 k )) f (-(A n + AA(0i k ))x + A m x + B m r). (5) 
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To take into account coupling between longitudinal and lateral dynamics due to the loss of symmetry in the case of 
damage, possible coupling nonlinearities are introduced into the model of F/A-18 dynamics. The resulting additional 
terms are combined into a single vector, namely as 


x — (A J; + AA($i))x + ( B n + A_8(6b))M + £ 


(6) 


where 


f = 


0 

-9a - 22.5/3 - 4.5 p -9 q— 18r - 4.5 pq - 4.5 pr - 9 qr 
0 
0 
0 
0 

-22.5 a - 13.5 /3 - 4.5 p - 13. 5g - 36r - 9 pq - 9 pr - A.5qr 
-4.5a - 22.5/3 - 13.5/7 - 22.5 q - 9 r - 18 pq - 18 pr - 4.5 qr 

0 

0 


(7) 


The effect of £ can be addressed using variable structure or sliding mode control for the compensation of the effect 
of nonlinearities. To accomplish this, ^ is included in the current model, with k corresponding to the current model, 
such that 

Xi = (A„ + AA0ik))x + (B m + AB02k))u + - Ae^. (8) 

Note that the vector estimating the nonlinearities is adjusted using 


h = -ei 


where E = diag[(sign(e^))] and 


1 = 


0 

22.5(|a| + ^61 + \p\ + \q\ + \r\) 2 
0 
0 
0 
0 


36(|ar| + \/3\ + \ P \ + \q\ + M) 2 
22.5(|a| + ^81 + \p\ + \q\ + |r|) 2 
0 
0 

The estimate of the nonlinearities is then added to the control law to give 

u = (B„ + A B02k))\-(A n + A A(6 lk ))x - + A m x + B m r). 


(9) 


( 10 ) 


As a next step in our evaluation of the suitability of the MMST approach to handle structural damage, we integrated 
the F/A-18 damage model with that of the TAMU wing, 20 and assumed coupling between these subsystems as well as 
parameter changes due to damage. The criterion for choosing the coupling terms was that the system quickly becomes 
unstable in the case of damage. 

The augmented model is of the form: 


xs - A_y (6\ )x s + Bs (62)11 + A R s(0\ )xr (11) 

Xr = Ar{0\)Xr + Br(02)u + Asr(0i)Xs + 1 (12) 

where x* e R" is the vector of rigid body states, x,y e R ; is the vector of pitch-and-plunge states, and u e R'" is the 
control input vector. As , Bs,Ars , Ar, B r , and As r are matrices of appropriate dimensions. It is assumed that the above 
system is completely controllable and observable, which is a realistic assumption. 

4 of 25 


American Institute of Aeronautics and Astronautics 



The control objective is to design a control input u(t) so that xr(J) follows a rigid-body reference model, while the 
pitch-and-plunge states are stabilized around zero, i.e. lim, xs (f) = 0 despite the large changes in 9 and £ : . 

We note that the baseline controller for this case is similar to the previous one, except that it is now designed for an 
extended state space, consisting of both rigid-body and structural states. A solution to this problem in the case of large 
uncertainties in 9 and xi is proposed in this paper, and consists of model set design and reduction, and implementation 
of a MMST scheme for the augmented damage model. 

III. A Model Set Design Approach 

In this section the focus is on the model set design for multiple model control. 

Let a parametric set be defined as: Sg = {6 : (0, ■)„„■„ < 9 , < ( 9 j) max , for i = 1,2 , ...,1V} where N is the number 
of uncertain parameters. Then the multiple model control technique involves placing a finite number of models within 
this set, designing corresponding observers, and on-line selection of a single controller from a set of controllers also 
corresponding to the predetermined models and observers. This selection is usually based upon some comparison 
metric to determine which observer is closest in some sense to the true system from observations made. A control 
strategy is then applied given the current observer selection. 

The choice of the control strategy is relevant for model set design based on performance specifications for a closed- 
loop system. This will be discussed in more detail below. 
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Figure 2. Example model set design with equally spaced models and the same number of model divisions for each uncertain parameter. 

For the F/A-18 application, it is assumed that the uncertain set is defined as Sg = {(ffi, 82) : -1 < 9 \ < 1; -1 < 
82 < 0}, where a single model can be specified by a point 8 = [ 8\ , 82 ] r € Sg. One possibility for the model set design 
is based on equally spaced models with the same number of model divisions for each uncertain parameter. An example 
of this approach is shown in Figure 2. For control the inverse dynamics strategy is chosen for this application, yielding 

u = (B„ + AB(fi 2 k)Y l (-(A n + A A(8 xk ))x + A m x + B„,r ) (13) 

where k is the model that is currently selected, and 8 \k and 82 k are the corresponding parameters of that model. This 
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strategy results in the following closed-loop system: 

x = [A„ + A A(0j) + (B„ + AB(0 2 ))(B n + AB(6 2k ))~ ] (-(A n + A A0 lk )) + A m )]x 

+(B„ + A B(02k)KB n + A B(0 2k ))- 1 B m r. (14) 

Proposed Model Set Design Approach: Our model set design approach is based on equally spaced models within 
the parametric set. We first define the “centered” models. For instance, if the model set has a single parameter such 
that 8 6 [0, 1], let TV be the number of models that will be used to cover this interval. Then the length of each 
subinterval will be 1 /A. Let us assume that model centers should be also placed at interval boundaries. In that case 
the total number of models is A + 1 , and the location of the “centered” models will be at 0, 1 /A, 2/A, ... 1-1 /A, 1 . 
If the models do not need to be placed at the interval boundaries, then there will be exactly A models located at: 
l/A/2, 3/A/2, ... 1 - l/A/2. 

The next issue to be addressed is how to choose A. One possibility is to make that choice based on some stability 
and/or performance criterion. We note that a given A defines a set of equally spaced centered models. For the 2 
parameter application it can be seen that for model j, the set of values that satisfy 

\\0-0jh < \\e-8ih for i + j (15) 

defines Sg p a rectangular subset of Sg for which model j is the closest. Observe that the controller for model j is 
defined by (13) where k = j yielding the closed-loop system in (14). To evaluate the case when j is the selected model 
and 8 e Sg. we consider the cases when the true plant is at the corners of Sg j . and evaluate if the resulting closed-loop 
systems satisfy performance specifications. For instance, one possible criterion is that the closed-loop eigenvalues 
must be within a specified euclidean distance from the desired eigenvalues of the reference model. By evaluating the 
corners of the models, one can ensure that, if the correct model is selected and the controller based on the model center 
is chosen, the closed-loop system will be guaranteed to be close to the reference model in the sense of eigenvalue 
location. Figure 3 illustrates this concept. 



A simple algorithm for ensuring the least number of models where all of the models meet the specifications de- 
scribed above is as follows: 

Algorithm: 

1 . Set maximum number of model divisions for each parameter 

2. Divide the space into even models 

3. Evaluate the closed-loop system with the control parameters at the model centers and the true system parameters 
at the model corners 

• If all models pass, then reduce the number of model divisions for each parameter and return to step 2 

• If any model fails, then increase the number of model divisions for each parameter to the value that passed 
in the previous iteration, and quit. 
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IV. Multiple-Model-based Control Design 


In this section we first give a general description of the Multiple Models, Switching and Tuning (MMST) technique, 
followed by a description of several different switching strategies that improve the overall system robustness. 

A. Multiple Models, Switching and Tuning (MMST) 

The concept of MMST is based on the idea of describing the dynamics of the system using different dynamic models for 
different operating regimes; such models are close in some sense to the dynamics of the system in different operating 
regimes, and the corresponding observers are consequently referred to as the identification models. The basic idea is to 
set up such identification models and corresponding controllers in parallel, Figure 4, and to devise a suitable strategy 
for switching among the controllers to achieve the desired control objective. While the plant is being controlled using 
one of these controllers, the identification models are run in parallel to generate some measure of the corresponding 
identification errors and to ultimately determin which observer is, in some sense, closest to the current operating regime 
of the plant. Once such an observer is found, the switching mechanism switches to (or stays at) the corresponding 
controller, where the switching interval is a parameter chosen by the designer. 

The main feature of this approach is that in LTI systems, if one of the models is adaptive where estimated parameters 
are limited only by the bounds of the uncertainty space, the approach results in a stable overall system in which 
asymptotic convergence of the output error to zero is guaranteed under relatively mild conditions. In addition, if 
adaptive control is used and the controller parameters are adjusted using adaptive algorithms with projection, the 
overall system is robust to different perturbations, including bounded external disturbances, time variations of plant 
parameters, and some classes of unmodeled dynamics. As also shown through extensive simulations, the performance 
of the overall switching system can be dramatically improved as compared to that achieved using a single adaptive 
controller. 



ADAPTIVE 

ESTIMATOR 

BANK 


CONTROLLER 

BANK 


Figure 4. Structure of the Multiple Model-Based Controller: Outputs of the parallel observers 0\, O 2 , ... 0,y are used to find that closest 
in some sense to the current plant dynamics, and switch to the corresponding controller. 

In the context of reconfigurable control design in the presence of parametric uncertainties and/or sensor, actuator 
and structural failures, the identification models (observers) 0 \ , ..., On from Figure 4 correspond to different regions in 
the parameter space characterizing different types of failures, while C 1 , ..., C’,v denote the corresponding controllers. 

In the context of damage accommodation, we use the MMST approach with a model that corresponds to the 
nominal (no-damage) case, and damage models that cover the entire parametric set. The best estimate at every instant, 
as selected by the MMST subsystem, is used in the control law. 

B. Switching Strategies 

Switching from one controller to another within the MMST framework is based on the comparison of estimation errors 
representing a difference between the outputs of on-line observers and observations of the the true system. For the 
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( 16 ) 


F/A-18 model the observers differ by the values of 0\ t and 62 /, and are of the form: 

Xj = (A„ + AA{6y))x + (B m + AB02i))u ~ d<?, for i - 1,2, ...n 


where n is the number of models, and e, = (x, - x). 

The metric employed for model comparison is the integral of the squared observer error over a finite time period. 
Corresponding performance indices are of the form: 


hit) = 



ej (T)e,(T)dT 


for i = 1,2, ...n. 


(17) 


and the current model selection index satisfies 


i* = arg min /,(f), for i = 1, 2, ...n. (18) 

i 

One disadvantage of the decision making based on this criterion is that, in the case of all fixed observers where 
all model parameter estimates are fixed at the model centers, switching to a wrong observer can occur. To gain more 
insight into this problem, the error model using equations (1) and (16), is written as: 


ei = Cl<pi - Ae, ( 19 ) 

where 0/ = [By - By 0 2 i - 62V and O = [ a 2 q bfi LE F ] yielding 

e,- = O.F<pi + & (20) 

Cly = — /tO/r + 0 (21) 

ki = -At,-n F fa ( 22 ) 


where if </>,■ is constant £■', - 0, while if 0/ is not constant but changes only infrequently <(',■ is an exponentially decaying 
signal whose decay rate is determined by the observer gain A. Now assuming that a 0 the performance indices are 
rewritten as: f 

l,(t) = f (T)e,(T)dT = <p[ ( f Qp(T)Q F (T)dTl^i. (23) 

Jt-At \Jt-At ! 

From this point one can see that in order to guarantee that the system will switch to model closest to the current plant 
dynamics in the euclidean sense (i.e. i* = arg min || (f) l H 2 ), the singular values of the matrix 


f 

Jt- A 


Qi(r)Qf(T)dr 


must be equal. Geometrically in the 2 parameter application the level sets of 


f 

Jr- A 


ej ( T)ei(r)dT 


(24) 


in the two dimensional parameter space defined by Sg must be circular in order to guarantee that model with the 
smallest parametric error will be selected. From (23) it can be seen that this coinsides with the singular values of (24) 
being equal. Since this is not guaranteed, switching to the wrong model can occur. From simulations and testing it 
has been found that, to reduce this undesired phenomena in multiple model selection, several techniques can be applied. 


Making the current model adaptive and all other models fixed: This method of all tested appears to have the great- 
est influence on undesired switching. For the current model estimates, the parameters are adjusted using the adaptation 
law 

0, = Proj^ | (25) 

where F = diag([yi , J2]), 71,72 > 0. Also note that Proj^ [(•)] is the projection operator with Sg j = {By, (b/) : By < 
B\i < By; 82 i < 62 i < By], where 0 ; , and 0 ;l - for j = 1,2 denote the lower and upper parameter boundaries respectively 
for an individual model i. This method reduces incorrect switching as follows. Through adaptation, the observer that is 
currently selected is able to outperform the other fixed observers by locating local parameters within the current model 
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boundaries. This typically occurs until the observer of a model that is truly closer to the actual system (i.e. in the sense 
of parameter location) results in a better performance index. 


Bias the current model : In this case the performance indices are modified as follows: 

?' (T)ej(T)dT -6 if i - k 


m = 


f i 

f H 


(T)ei(r)dT 


otherwise 


(26) 


where 6 > 0 is the bias, /, is the comparison metric, and k is the current model. With this technique the selection will 
favor the current model unless the current model error is sufficiently large. 


Decrease the selection frequency with adaptation of the current model. Selection frequency has also been found to re- 
duce unnecessary model switching when the current model is adaptive and other models are fixed. By increasing the 
time between model selection it is assumed that the adaptive model has more opportunity to adapt and outperform other 
incorrect fixed models. 


V. Simulations 

Simulations of the multiple model design and control strategies described is given below. 

The simulation scenario is the F/A-18 performing a pitch doublet maneuver. It is assumed that the damage pa- 
rameters change from 8 = [0, 0] r to 8 = [0.4687, - 0.4653 ] T at 3 seconds. For model set design, the closed loop 
eigenvalues are constrained to be within 0.2 of the reference model eigenvalues. This leads to a set of equally spaced 
models with five divisions for each parameter, or the total of 25 models. First the results for all fixed models with an 
observer gain of A = 20 and integration time of At = 0.2 seconds are given. 

A. Damage model without augmentation 

This case served primarily for choosing the most suitable MMST strategy and tuning of the MMST controller. The 
response of the system with the MMST controller is shown in Figures 5 and 6. It can be seen that with fixed models 
the system states are well behaved but incorrect switching occurs. 

This however, is shown to be remedied in the following simulation where switching improvement techniques de- 
scribed earlier are employed. Under the same scenario, the current model is made adaptive, while all others are fixed. 
Adaptation gains are set to j\ = 70 and ys = 500, the observer gain is set to A = 20, and At = 0.2. Also, biasing is 
applied with 

6(t) — 0.2 I ej (r)e,'(T)c/T 
Jt-At 

and the switching frequency is reduced to wait 0.2 seconds between model selections. Results for this simulation are 
given in Figures 7 and 8. The longitudinal states are shown to track the reference model accurately. It is also seen that 
the switching is well behaved, due to the use of the previously described techniques that prevent switching to a wrong 
model. 

B. Damage model with inter-axes coupling 

Without compensation, the addition of these nonlinear terms to the F/A-18 dynamics results in instability as shown 
in Figure 9. In this simulation the same MMST controller described above (with the current model adaptive, biasing, 
and decreased switching frequency) is applied to the damage scenario where, after 3 seconds, 8 = [0, 0 ] 7 changes to 
8 = [0.4687, - 0.4653] r , and f is introduced. 

With the addition of the variable structure term to the previous multiple model controller, the same highly nonlinear 
system can be stabilized, and can accurately track the reference system. This is shown in Figure 10, which shows the 
longitudinal states of the nonlinear system throughout the doublet command. Figure 12 shows the actual nonlinearities 
and the variable structure terms used for cancellation. The damage parameters of the system can also be seen in 
Figure 1 1 where the model switching and parameter estimates are well behaved with the addition of fk to the observers 
and controllers. 
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time [s] 


Figure 5. MMST with Fixed Models: Response of the longitudinal states under damage at t = 3 seconds. 


C. MMST-based Control of an Aero-Servo-Elastic Model 

We first simulated the case when the system is controlled using a reduced-order controller, i.e. the controller that 
does not take into account pitch-and-plunge states. Figure 13- 14(b) show the resulting response. It is seen that 
the longitudinal state errors are large, while the pitch-and-plunge states are highly oscillatory. The response of the 
parameter estimates is also seen to be highly oscillatory, while the switching pattern exhibits a similar response. 

The next simulation is of the case when the system is controlled using a full-order controller, i.e. the controller 
that takes into account pitch-and-plunge states. Figure 1 6—1 7(b) show the resulting response. It is seen that the 
longitudinal state errors are much smaller than in the previous case, while the oscillations of the pitch-and-plunge 
states are suppressed. The system switches only few times appropriately despite damage, active structural modes, and 
coupling nonlinearities. 

VI. MMST-based Guaranteed Performance Adaptive Control (GPAC) 

One of the fundamental questions regarding the proposed damage-adaptive control scheme is that of guarantees 
that can be given that the proposed system will behave in the desired way under all possible values of 6 and xi. 

This is consistent with one of the main challenges in adaptive control, where an open issue is to derive suitable 
performance metrics for adaptive system, and develop strategies to assure that the dynamics of the overall adaptive 
system is close to that of a reference model for all time. If such performance metrics and control strategies could be 
found, then one could guarantee that both the transient and steady-state properties of the closed-loop system are close 
to that of the reference model. The most challenging issue here is that this needs to be achieved for any variation in 
plant parameters. 

Here we present an initial analysis of the problem of achieving consistent system performance in adaptive control, 
and propose several solutions based on MMST. 

Let the plant dynamics be described by: 


x = f(p, x, u), 


(27) 
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(a) 




(b) 


Figure 6. MMST with Fixed Models: Response of the estimates of the damage parameters under damage at / = 3 seconds. 
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Figure 7. MMST with one Adaptive Model, Biasing and Decreased Selection Frequency: Response of the longitudinal states under damage 
at t = 3 seconds. 


where x e R — » K", u e R — » IR”', and p e S p , where S p is a known set. 

Reference Model: Let a stable reference model be given as: A* = g(x*, r), where r is a bounded piece-wise 
continuous reference input. 

Adaptive Controller: Let u(t) be an output of the adaptive controller defined as: 

u = h(x, r, p), (28) 

where p denotes an estimate of p generated by: 

x = f(p, x, u ) + A(x - x), p = q(x, u, x), (29) 


and let u* = /z(x*, r, p). 

We now consider the following definitions: 

Definition 1 (Consistent Instantaneous Performance): The closed-loop system (27)-(29) achieves consistent instan- 
taneous performance if and only if, for all ||.r 0 || < 6 X , all ||x*|| < 6 X , all p e S pi and given 5 X > 0 and 6 U > 0, 

\m - **(011 < s x , \\u(t) - u*(t)\\ < 6 U , Vf > t 0 . (30) 

Definition 2 (Consistent Interval Performance): The closed-loop system (27)-(29) achieves consistent performance 
over an interval [t„, T ] if and only if, for all x„ and x* in some predefined sets of initial conditions, all p e S p . and 
given e x and e„. 



l|x(0 - x\t)\\dt < e x . 



|| u(t) - u(t)\\dt < 


e u - 


(31) 


The main challenge here is to design a control strategy that will keep the above measures at or below their pre- 
specified bounds for any change in the plant parameters. This is a highly complex problem whose solution depends on 
many factors. Some possible solutions are proposed below. 
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Figure 8. MMST with one Adaptive Model, Biasing and Decreased Selection Frequency: Response of the damage parameter estimates 
under damage at t = 3 seconds. 
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Figure 9. Response of the longitudinal states in the case of damage with coupling terms included, and no compensation. 
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Figure 10. Response of the longitudinal states with coupling terms included and variable structure terms used to compensate for nonlinear- 
ities. 
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(a) 




(b) 

Figure 11. Response of the damage parameters and their estimates with variable structure terms used to compensate for nonlinearities. 
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Figure 12. Response of the coupling terms and their estimates with variable structure terms used to compensate for nonlinearities. 



Figure 13. State response of the ASE model to damage at t = 3 seconds in the case without compensation. 
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(a) 




(b) 


Figure 14. Response of the damage parameters and their estimates with no compensation for structural dynamics. 
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Figure 15. Response of the pitch-and-plunge states to damage at t = 3 seconds in the case without compensation. 



Figure 16. State response of the ASE model to damage at t = 3 seconds in the case with compensation. 
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Figure 17. Response of the damage parameters and their estimates with compensation for structural dynamics. 
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Figure 18. Response of the pitch-and-plunge structural states to damage at t = 3 seconds in the case with compensation. 


MMST with Model Injection: In this case the idea is to use all fixed observers and populate the subset to which the 
system switches with a large number of observers. Figure 19. This way at least one of the observers will be fairly 
close to the true plant dynamics. If the system adapts from there, the tracking error can be minimized along with the 
Consistent Performance performance index. 

Simulation of a first-order plant with an uncertain control derivative coefficient is used to evaluate the model- 
injection technique, and the resulting values of the Consistent Performance Index (CPI) against the level of uncertainty 
(from 0.05 to 1) is shown in Figure 20. It is seen that the model injection technique has a potential to minimize the CPI 
for any value of the uncertain parameter. 

MMST with Adaptive Gain Scheduling: In this case all the observers are adaptive but the parameter estimates are 
adjusted using different adaptive gains, and the system will switch to the observer with the adaptive gain that is ’’the 
best” for the current operating regime. This approach is illustrated in Figure 21. 

The same first-order example is used to illustrate the adaptive gain scheduling approach, and the results are shown 
in Figure 22. It is seen that the values of CPI are similar to those obtained with model injection. Hence this is another 
potentially feasible approach to CPI minimization. 

MMST design using Model Sets: This strategy is based on the use of an output error feedback (OEF) term in the 
control law, Figure 23, and robust stability results to stabilize all plants within a subset of the parametric set. As shown 
recently, 1 the resulting closed-loop adaptive control system is linear and time-varying. This OEF design procedure can 
be used for the Model Set Design, which addresses the question as to how many models should be used, and where they 
should be placed within the parametric set. The proposed approach is based on extending the current set until the robust 
stability conditions are violated, and designing a controller with output error feedback for the next subset. An important 
aspect of this approach is that, when the system switches to the right set, the resulting closed-loop system will have 
some minimum performance specified during the model-set design. For this reason this approach is also referred to as 
the Guaranteed Performance Models Set Design (GP-MSD). Simulation results using this strategy are shown in Figure 
24 resulting in performance similar to that obtained in the previous cases. These strategies will be compared to find the 
one most suitable for minimizing the CPI in the case of multivariable plants with multiple uncertain parameters. 
Guaranteed switching strategies: One of the major questions in MMST design is that of convergence of the switching 
parameter to the value closest to or corresponding to the actual plant dynamics. Many theoretical issues need to be 
addressed before one can derive analytic conditions for a finite-time interval over which the system is guaranteed 
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Figure 19. Strategy with model injection: The system will first switch to a fixed observer M 4 closest to the current plant dynamics. At that instant 
the set corresponding to M 4 is is populated with multiple observers. Adaptation is hence replaced by rapid switching to the observer closest to the 
current plant dynamics. 


to switch to the ’’right” model. These include: performance criterion used for switching; model set design; overall 
switching strategy; adaptive control law; speed of parameter variations in the system; and signal excitation. 

This aspect of the MMST design is extremely important since switching to the right model in finite time guarantees 
some minimum performance. Theoretical and practical aspects of this strategy need to be investigated in detail. 

VII. Conclusions 

In this paper we present a comprehensive approach to damage accommodation using Model Set Design, MMST, 
and Variable Structure compensation for coupling nonlinearities. 

In the paper a new approach is described related to the design of model sets for the MMST, given the bounds on the 
plant parameters. The MMST controller also includes terms that compensate for cross coupling nonlinearities arising 
due to the damage. 

We also implemented the MMST approach to an integrated rigid-body and structural dynamics model that combines 
linearized F/A- 1 8 dynamics with the pitch-and-plunge dynamics of the TAMU wing. 

In all cases excellent performance was obtained despite damage, coupling nonlinearities, and structural dynamics. 
This demonstrates the feasibility of the MMST approach for damage accommodation. 

While the augmented ASE model used in this paper is simple and based on ad-hoc assumptions, our further work 
will include testing of the IDACS controller on more realistic models such as the Generic Transport Model (GTM) 
simulation that accounts for changes in aerodynamics due to damage, as well as for changes in the center of gravity, 
retrofitted with an aeroservoelastic model. 
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Figure 21. Different adaptive gains will result in different system trajectories (left figure). The MMST approach will result in moving along a 
current trajectory ( solid blue line ) or switching from one trajectory to another ( dashed blue line ) to improve the response ( right figure). 
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Figure 22. Comparison of the strategy with adaptive gain scheduling with the following strategies: all-fixed observers, all-fixed plus one adaptive 
observer, and all-fixed plus one adaptive observer with scheduled gains. 



Figure 23. Structure of the system with output error feedback. 
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Figure 24. Comparison of the GP-MSD strategies. 
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